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I. INTRODUCTION 

The integer quantum Hall effect in two-dimensional 
(2D) electron gas arises from the quantized 2D Landau 
levels (LL). The non-trivial band structure topology is 
characterized by non-zero Chern numbers 2 ^ 3 -. Later on, 
quantum anomalous Hall insulators based on Bloch-wave 
band structures were proposed in the absence of Landau 
levels^. In recent years, the study of topological insula- 
tors (TI) in both 2D and three dimensions (3D) has be- 
come a major focus of condensed matter physics^—. TIs 
maintain time-reversal (TR) symmetry, and their band 
structures are characterized by the nontrivial .^-indcx. 
As for the 3D TIs, various materials with Bloch-wave 
band structures have been realized and the stable he- 
lical surface modes have been detected^— . Since LL 
wavefunctions have explicit forms with elegant analyti- 
cal properties, TIs based on high dimensional LL struc- 
tures would provide a nice platform for further theoret- 
ical studies. In particular, interaction effects in the flat 
LLs are non-perturbative, which could lead to non-trivial 
many-body states in high dimensions. 

The seminal work by Zhang and Hu^i generalizes LLs 
to the compact S^-sphere with particles coupled to the 
SU(2) gauge potential. The isospin of particles I scales 
as R 2 where R is the radius of the sphere. Such a system 
realizes the four dimensional integer and fractional TIs. 
The 3D and 2D TIs can be constructed from the 4D TIs 
by dimensional reduction—,. Further generalizations to 
other manifold have also been developed^— . Recently, 
the LLs of non-relativistic fermions have been general- 
ized to arbitrary dimensional flat space R D by two of 
the authors^. For the simplest case of 3D, the SU(2) 
Ab.aranov-Cash.er gauge potential replaces the role of the 
usual U(l) vector potential. Depending on the sign of 
the coupling constant, the flat LLs are characterized by 
either positive or negative helicity. In the positive and 
negative helicity channels, the eigenvalues of spin-orbit 
coupling term a ■ L take values of I and — (I + 1), respec- 
tively. Each LL contributes a branch of helical surface 
modes at the open boundary. When odd numbers of LLs 



are fully filled, there are odd numbers of helical Fermi 
surfaces. Thus the system is a 3D strong topological 
insulator. This construction can be easily generalized 
to arbitrary D-dimensions by coupling the fundamental 
spinors to the SO(D) gauge potential. 

Quantized LLs of 2D Dirac fermions have also been ex- 
tensively investigated in the field theory context known 
as the parity anomaly^Sr— . This can be viewed as the 
square root problem of the usual 2D non-relativistic LLs. 
External magnetic fields induce vacuum charges with the 
density proportional to the field strength. The sign of the 
charge density is related to the sign of the fermion mass. 
There is an ambiguity if the Dirac fermions are massless. 
In this case, there appear a branch of zero energy Lan- 
dau levels. Each of them contributes ±g fermion charge. 
It is similar to the soliton charge in the Jackiw-Rebbi 
mode l 34 i 35 , which is realized in condensed matter systems 
of one dimensional conducting polymer 36 . Depending on 
whether the zero energy Landau levels are fully occupied 
or empty, the vacuum charge density is wnere V 

is the magnetic length. In condensed matter physics, the 
best known example of Dirac fermions is in graphene, 
which realizes a pair of Dirac cones. The quantized LLs in 
graphene have been observed which distribute symmetri- 
cally with respect to zero energy. Their energies scale as 
the square root of the Landau level index. The observed 
Hall conductance per spin component are quantized at 
odd integer values, which reflects the nature of two Dirac 
cones in graphene for each spin component^—. 

In this article, we generalize the LLs with full rota- 
tional symmetry of Dirac fermions to the three dimen- 
sional flat space and above. It is a square root problem 
of the high dimensional LLs investigated in Ref. [l| . Our 
Hamiltonian is very simple: replacing the momentum op- 
erator in the Dirac equation by the creation or annihila- 
tion phonon operators, which are complex combinations 
of momenta and coordinates. The LLs exhibit the same 
spectra as those in the 2D case but with the full rota- 
tional symmetry in D-dimensional space. Again the zero 
energy Landau levels are half-fermion modes. Each LL 
contributes to a branch of helical surface mode at open 
boundaries. 
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This paper is organized as follows. In Sect. [H] af- 
ter a brief review of the 2D LL Hamiltonian of Dirac 
fermions in graphene, we construct the 3D LL Hamilto- 
nian of Dirac fermions. Reducing this 3D system to 2D, it 
gives rise to 2D quantum spin Hall Hamiltonian of Dirac 
fermions with LLs. In Sect. HID we further solve this 3D 
LL Hamiltonian of Dirac fermions, and its edge proper- 
ties are discussed. For the later discussion of generalizing 
the 3D LL Hamiltonian to arbitrary higher dimensions, 
in Sect. IIV[ we briefly review some properties of D di- 
mensional spherical harmonics and spinors. In Sect. |V] 
and Sect. IVI1 we extend the solutions of LL Hamiltoni- 
ans to arbitrary odd and even dimensions, respectively. 
Conclusions are given in Sect. IVIII 



II. THE LANDAU LEVEL HAMILTONIAN OF 
3D DIRAC FERMIONS 

A. A Brief Review of the 2D LL Hamiltonian 

Before discussing the LL problem of Dirac fermions in 
3D, we briefly review the familiar 2D cas o 39 ' 40 to gain the 
insight on how to generalize it to high dimensions. The 
celebrated condensed matter system to realize 2D Dirac 
fermion is the monolayer of graphen o 39 ' 40 , which pos- 
sesses a pair of Dirac cones with spin degeneracy. Here 
for simplicity, we only consider a single 2D Dirac cone 
under a uniform magnetic field Bz. The Landau level 
Hamiltonian in the xy-plane reads 



H 



2D.LL 



vf^{Px - ^A x )a x + (p y - ^A y )a y }, (1) 



where the Dirac fermion with momentum p is minimally 
coupled to the U(l) magnetic field with symmetric gauge 
potentials 



I _ " A - B 



D 
~2 



(2) 



satisfying V x A = Bz; the Fermi velocity vp is related 
to the cyclotron frequency oj via the magnetic length /' 
as 




I'oj 



(3) 



For later convenience, we define Iq — y/21' which will be 
termed as cyclotron length below. The spectra of Eq. Q] 
consist of a branch of zero energy LL, and other LLs with 
positive and negative energies distribute symmetrically 
around zero energy. The energy of each LL scales as 
the square root of the Landau level index. It is well- 
known that Eq. [T]can be recast in term of creation and 
annihilation operators 



*12D,LL — —F= 
V2 
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it 



ia x 



(4) 



where di{i = x,y) are the phonon annihilation operators 
along the x and y-directions, with the form as 



a, 



1/1 ,.lo 



(5) 



In Eq. HI two sets of creation and annihilation oper- 
ators combine with 1 and the imaginary unit i. In or- 
der to generalize to 3D, in which there exist three sets 
of creation and annihilation operators, we employ Pauli 
matrices to match them as explained below. 



B. The construction of the 3D LL Hamiltonian 



We define the rotationally invariant operator B as 



B. 



3D 



J_ (Pik _ .n 
s/2 1 h 



(0) 



where the repeated index i runs over x, y and z; &i is the 
phonon annihilation operator along the i-direction; is 
the cyclotron length. We design the 3D Landau level 
Hamiltonian of Dirac fermions as 



H 3D = 



B 3D 



(7) 



Eq. [7] contains the complex combination of momenta 
and coordinates, thus it can be viewed as the generalized 
Dirac equation defined in the phase space. Using the 
convention of a, j3 and 7-matrices defined as 
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Eq. [7] is represented as 

H 3 d = vf ^ { a iPi + 7;*frp-} 



(8) 



i—x,y,z 



where vf = A mass term can be added into Eq. [S] 

as 



H. 



3D,ms 



= A/3 = 



AJ2X2 






-A/2X2 



(9) 



A similar Hamiltonian was studied before under the name 
of Dirac oscillato r 41 ' 42 , which corresponds to Eq. [7] 
plus t he mass term of Eq. |9] with the special relation 
lo = ^/hc 2 /Auj. However, the relation between the solu- 
tion of such a Hamiltonian to the LLs and its topological 
properties were not noticed before. 

The corresponding Lagrangian of Eq. [S] reads 



£ 



ip^-yoihdt - iv^hdi^ip - v F hipi^ ^iipF 0l (r), (10) 



3 



where F° l — Xi/l 2 ,. Compared with the usual way that 
Dirac fermions minimally couple to the U(l) gauge field, 
here they couple to the background field in Eq. [TUl 
through ijoji- It can be viewed as a type of non- minimal 
coupling, the Pauli coupling. Apparently, Eq. [7] is rota- 
tionally invariant. It is also time-reversal invariant, and 
the time-reversal operation T is defined as 



(11) 



where K represents the complex conjugation operation, 
and T 2 = -1. 



C. Reduction to the 2D quantum spin Hall 
Hamiltonian of Dirac fermions with LLs 

If we suppress the z-component part in the definition 
of Eq. [6j we will arrive at double copies of the usual LL 
problem of 2D Dirac fermions with Kramer degeneracy, 
which can be considered as the ^-topological insulator 
Hamiltonian arising from LLs of 2D Dirac fermions. We 
define the operator Bid as B2D — —i&xO<x — icyO-y, and 
the Eq. [7] reduces to 
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4 + ih l \ 



Vp 



it + ial 
jj < 1 r 
V a y — ia x J 

( p- - A- 

P++A+ 

p-+A_ 
V p+ - A+ 



V 



where p± — p x ± ip y and A± = A x ± iA y . It is reducible 
into a pair of 2 x 2 matrices as 



#21), ± — Vp 



p- ± A- 
P+ ±A+ 



(12) 



which are time-reversal partner to each other. Thus Eq. 
[T^] can be viewed as the quantum spin Hall Hamiltonian 
of 2D Dirac fermions. 

A similar situation occurs in the strained graphene sys- 
tems in which lattice distortions behave like a gauge field 
coupling. Signatures of LLs due to strains have been 
observed in Ref. |43j. Due to the TR symmetry, the 
Dirac cones at two non-equivalent vertices of the Bril- 
louin zone see gauge fields with a opposite sign to each 
other. Such a coupling is also spin- independent. How- 
ever, the TR transformation connecting two Dirac cones 
satisfies T 2 = 1, thus LLs due to strain are not topo- 
logically protected. They are unstable under inter-valley 
scattering. 

Eq. rT2] exhibits the standard minimal coupling to 
the background U(l) gauge field. Its solutions are well- 
known thus will not be repeated here. After all, there is 



no non-minimal coupling in 2D. Each state of the zero en- 
ergy LL is actually a half-fermion zero mode. Whether it 
is filled or empty contributes the fermion charge ± i . As 
the chemical potential ^ = , magnetic field pumps vac- 
uum charge density p(r) — ±|^-S. In the field theory 
context, this is an example of the parity anomaly^Sr— . 
Our 3D version and generalizations to arbitrary dimen- 
sions exhibit similar effects as will be discussed below. 



III. THE BULK SPECTRA OF THE 3D DIRAC 
FERMION LLS 

In this section, we will present the solution of the 
spectra and wavefunctions of the 3D LL Hamiltonian for 
Dirac fermions, which can be obtained basing on the so- 
lutions of the 3D LL problem of the non-relativistic case 1 . 
We start with a brief review of the non-relativistic case 
of the 3D LL problem. 



The 3D isotropic non-relativistic LL 
wavefunctions 



The 3D isotropic LL Hamiltonians for non-relativistic 
particles are just spin-i fermions in the 3D harmonic 
oscillator plus spin-orbit coupling^ as 



(13) 



Their eigenfunctions are essentially the same as those of 
the 3D harmonic oscillator of spin-^ fermions organized 
in the total angular momentum eigen-basis of j, j z as 



(14) 



where n r is the radial quantum number; j± — i± ^ repre- 
sent positive and negative helicity channels, respectively; 
and I is the orbital angular momentum. Please note that 
I is not an independent variable from j±. We write it 
explicitly in order to keep track of the orbital angular 
momentum. The radial wavefunction can be represented 
through the confluent hypergeometric functions as 



R nrl (r) = N nrJ (^) l F(-n r ,l + ^,^) 
\L / Z l Q 



(15) 



where F is the standard first kind confluent hypergeo- 
metric function, 



F(-n r ,l+-,j 2 -) = x 



n 



^ T(-n r ) T(l + l + n) 
1 



T(n + 1) \l 2 



(16) 



When n r is a positive integer, the sum over n is cut off 
at n r . The normalization factor reads as 



N, 



In 



Tl r ,l 



/ 2r(z + 7i r + §) 

r(i + §) V r(n r + 1) ' 



(17) 
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The angular part of the wavefunction is the standard 
spin-orbit coupled spinor spherical harmonic function, 
which reads as 



Depending on the sign of the spin-orbit coupling in Eq. 
f]~51 one of the two branches of positive or negative helicity 
states are dispersionless with respect to j, and thus are 
dispersionless LLs. For H^u- , the positive helicity states 
become dispersionless LLs as 



Hz D -i>n r ,j+,l,j z = (2n r + ^)f^^nr,j+,l,jz ( 19 ) 

where n r serves as Landau level index. However, the neg- 
ative helicity states are dispersive whose eigen-equation 
reads 



HzD-i>n r ,j-,l,jz = ( 2 «r + 2/ + -)&J^ nrJ _ Aj - 2 . (20) 

Similarly, we have the following eigen-equations for the 
i?3_D i+ as 

3 

HzD,+^n r ,j + ,l,3z = i 2n r +2Z+ -)^Vr lr ,j + ,y 2! -(21) 

In this case, the negative helicity states become disper- 
sionless LLs with respect to j , while the positive ones are 
dispersive. 



B. 3D LL wavefunctions of Dirac fermions 

Now we are ready to present the spectra and the four- 
component eigenfunctions of Eq. [5] for the massless case. 
Its square is block-diagonal, and two blocks become the 
non-relativistic 3D Landau level Hamiltonians with op- 
posite signs of spin-orbit coupling studied in Ref. 



H 3D 



H- 
H+ 



P i 1 71 ,r 2 2 

— Mu r 

2M 2 



>{£•*+§»} 



i o 
o -/ 



(22) 



where M is defined through the relation Iq = \J h/(Moj). 

Its eigenfunctions can be represented in terms of non- 
relativistic Landau levels of Eq. [14] as presented in Sect. 
IIII Al Eq. [7] has a conserved quantity as 



K 



a-l + h 
-(a-l + h) 



(23) 



According to its eigenvalues, the eigenfunctions of Eq. [7] 
are classified as 



- 11 *±n r ,j,jz 



(24) 



respectively. ^± nr jj r is dispersionless with respect to 
j, while ^± Ur jj z is dispersive, respectively. The disper- 
sionless branch is solved as 



*i ■ (r) 



1 ^n r j + ,i, 3 Jr) 
,/2 [ ±iip nr -i tj _ t i +1Jz (r) 



with the energy 



E- 



(25) 



(26) 



Please note that the upper and lower two components 
of Eq. [25] possess different values of orbital angular mo- 
menta. They exhibit opposite helicities of j±, respec- 
tively. The zeroth Landau level (n r — 0) is special: only 
the first two components are non-zero. 

On the other hand, the wavefunctions of the dispersive 
branch read 



ii 



" \[2 



+fynr,j-,l+l,j„(f) 



with the spectra solved as 



(27) 



(28) 



These states are just discrete energy levels lying between 
two adjacent LLs. For simplicity, let us only consider 
the positive energy states. The degeneracy of these mid- 
gap states lying between the n-th and (n + l)-th Landau 
levels with n = n r + j + \ is finite, n(n + 1), due to 
finite combinations of n r and j. In particular, between 
the zeroth LL and the LLs with n r — ±1, these discrete 
states do not exist at all. 

Because Eq. [S] satisfies ^H^dP = —H$d, its spectra 
are symmetric with respect to zero energy. If the ze- 
roth branch of Landau levels (n r — 0) are occupied, each 
of them contribute a half-fermion charge. The vacuum 
charge is p 3D (r) = \ (f)*0y,j (j)> which are 

calculated as 



„3D 



(r) 



= Xt 



l + l 



r 2, I, 



2/3 I 2tt 



2=0 



T(l + I) u 



2irl 



2-kII' 



(*)K* 



as r 



(29) 



In 2D, the induced vacuum charge density in the gap- 
less Dirac LL problem is a constant, p 2D (r) = -^p = 

j^-^B , which is known as "parity anomaly". However, 
in the 3D case, the vacuum charge density p 3D (r) di- 
verges linearly, which is dramatically different from that 
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in 2D. This can be easily understood in the semi-classic 
picture. Each Landau level with orbital angular momen- 
tum I has a classic radius r; — V211q. In 2D, between 
T\ < r < rj-|_i, there is only one state. However, in 3D 
there is the 2j+ + 1 = 21 + 2 fold degeneracy, which is 
the origin of the divergence of the vacuum charge density 
as r approaches infinity. Generally speaking, in the case 
of D dimensions, the degeneracy density scales as r D ~ 2 
as shown in Sect. IVland Sect lVIl The intrinsic difference 
between high-Z? and 2D is that the high dimensional LL 
problems exhibit the form of non-minimal coupling. In 
2D, due to the specialty of Pauli matrices, this kind of 
coupling reduces back to the usual minimal coupling. In 
Eq. [TUl the background field is actually a linear divergent 
electric field, not the magnetic vector potential. Eq. [29] 
can be viewed as a generalization of "parity anomaly" to 
3D for non-minimal couplings. 

Now we consider the full Hamiltonian with the mass 
term Eq. [§] The mass term mixes the LLs in Eq. [55] 
with opposite level indices ±n r but the same values of j 
and j z . The new eigenfunctions become 



i'J 



-n T ,3,3* 



cost 
sin t 



— sin 9 
cos 9 



where cos 2 9 = \ 
tra are 



l + Vn r /Vnr + [A/(fiw)] 5 



, (30) 



The spec- 



E 



±\/n r (huj) 2 + A 2 . 



(31) 



The zeroth LL \?„ r _o j j (?) singles out, which is not af- 
fected by the mass term. Only its energy is shifted to 
A. 



C. Gapless surface modes 

As shown in Eq. (2H the 3D LL system of Dirac 
fermions has a center, and docs not have translational 
symmetry. Thus how to calculate its topological in- 
dex remains a challenging problem. Nevertheless, we 
can still demonstrate its non-trivial topological proper- 
ties through the solution of its gapless surface modes. 

We consider the surface spectra at a spherical bound- 
ary with a large radius R 3> ^o- The Hamiltonian H r< R 
inside the sphere takes the massless form of Eq. [8] while 
Hr>R outside takes the mass term of Eq. |H] in the limit 
of |A| —> oo. Again the square of this Hamiltonian 
(Hr<R+Hr>ii) 2 is just Eq. [52]subject to the open bound- 
ary condition at the radius of R. The spectra of the open 
surface problem of the non-relativistic 3D LL Hamilto- 
nian have been calculated and presented in Fig. 3 in Ref. 
[l]|, in which the spectra of each Landau level remain 
flat for bulk states and develop upturn dispersions as in- 
creasing j near the surface. The solution to the Dirac 
spectra is just to take the square root. Except the zeroth 
LL, each of the non-relativistic LL and its surface branch 



split into a pair of bulk and surface branches in the rel- 
ativistic case. The relativistic spectra take the positive 
and negative square roots of the non-relativistic spectra, 
respectively. The zeroth LL branch singles out. We can 
only take either the positive or negative square root, but 
not both. It surface spectra are upturn or downturn with 
respect to j depending on the sign of the vacuum mass. 
For the current Hamiltonian, only the first two compo- 
nents of the zeroth LL wavefunction are non-zero, thus 
it only senses the upper 2x2 diagonal block of vacuum 
mass in Eq. [5] thus its surface spectra are pushed upturn. 



IV. REVIEW OF D-DIMENSIONAL 
SPHERICAL HARMONICS AND SPINORS 

We will study the LL of Dirac fermions for general 
dimensions in the rest part of the paper. For later con- 
venience, we present here some background knowledge 
of the SO(D) group which can found in standard group 
theory textbooks^. 

The _D-dimensional spherical harmonic functions 
Yii m \(Cl) form the representation of the SO(D) group 
with the one-row Young pattern, where I is the number 
of boxes and {m} represents a set of D — 2 quantum 
numbers of the subgroup chain from SO(D — 1) down to 
SO (2). The degeneracy of Y/ : { m } is 



d m (SO(D)) = (D + 2l-2) 



(£>■ 



3)! 



l\(D-2)\ 



(32) 



Its Casimir is Y^i<j ^ij 



l(l + D-2)h 2 , where the orbital 
angular momenta are defined as = r^-pj — rjPi. 

We also need to employ the r-matrices. The 2x2 
Pauli matrices are just the rank-1 r-matrices. They can 
be generalized to rank-fc r-matrices which contains 2k + l 
matrices anti-commuting with each other. Their dimen- 
sions are 2 k x 2 k . A convenient recursive definition is 
constructed based on the rank- (A; — 1) r-matrices as 



r (fe) 

i 

2k+l 



r (k-l) 

1 a 



r (fc-i; 

J- a 

i o 

-I 



,(fc) 

2A- 



-il 
il 



(33) 



where i = 1, 2k— 1. For D = 2fc+l-dimensional space, 
its fundamental spinor is 2 fe -dimensional. The generators 



are constructed Si. 



,(fe) 



where 



(34) 



For the D = 2fc-dimensional space, there are two irre- 
ducible fundamental representations with 2 k ~ 1 compo- 
nents. Their generators are denoted as Sy and S'^, re- 
spectively, which can be constructed based on both rank- 

(fe-i) r 



(fe-i) 



and T 



(fe-i) 



matrices. For the first 2k — 1 



6 



dimensions, the generators share the same form as 



Sin — Si J 



1. 



-r 



(k-i) 



1 < i < j < 2k - 1, (35) 



"13 "tj 2 " *j 
while other generators Si^fc and S| 2 fe differ by a sign as 
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2fc 



±-rp 1] , l < i < 2fe — l. 

2 — — 



(36) 



We couple ^,{ m } to the SO{D) fundamental irre- 
ducible spinors. For simplicity, we use the same symbol s 
in this paragraph to denote the fundamental spinor repre- 
sentation (Rep.) for SO(D) with D = 2k + 1 and the two 
irreducible spinor Reps, for SO(D) with D = 2k. The 
states split into the positive (J+) and negative hc- 
licity sectors. The bases are expressed as yj ± ; S ,i.{j m }(&), 
where {j m } is a set of D—2 quantum numbers for the sub- 
group chain. The degeneracy number of D^ + ; S ,i,{,j m }(^) 
is 



d- =d {D + l ~ 2) - 
3+ s U(D-2)\ ' 



where d s is the dimension of the fundamental spinor 
representation. Similarly, the degeneracy number of 



d- -d ( D+l -w 

J - s {l-l)\{D-2)V 



(38) 



The eigenvalues of the spin-orbit coupling term 
V ( . ■ \ for the sectors of 3j + ; S ,i,{j m }(^) and 

yj-;s,l,{j m }(&) are I H and — (l + D — 2)fi, respectively. We 
present the eigenstates of the _D-dimensional harmonic 
oscillator with fundamental spinors in the total angular 
momentum basis as 

<. J± , s ,U, m }(r) = Rn r ,l(r)y j± ,s Mjm }(ti), (39) 
where the radial wavefunction reads 

RnM = K.i (£) e~^F(-n r , l + j,^) (40) 
and the normalization constant reads 

_D_ 
/ 2 



2T(n r + l + %) 



T(l + ID) V r(n r + 1) 



(41) 



V. THE LLS OF ODD DIMENSIONAL DIRAC 
FERMIONS 

In this section, we generalize the 3D LL Hamiltonian 
for Dirac fermions to an arbitrary odd spatial dimensions 
D = 2k + 1. We need to use the rank-fc T-matrices, 
which contains 2k + 1 anti-commutable matrices at the 
dimensions of 2 k x 2 k denoted as r| fc) (l < i < 2k + 1). 

(k) 

The definition of T,- and the background information of 



the representation of the SO(D) group is given in Sect. 

EH 

(k) 

We define B2k+i = <k and the 2k + 1- 

dimensional LL Hamiltonian of Dirac fermions i?2fe+i in 
the same way as in Eq. [71 Again the square of i?2fe+i 
reduces to a block-diagonal form as 



{H2k+lY 



+ 



p 

2M 
2fc4 



1 



+ -MwV - u 



(fc) 



In} 



I 
-I 



(42) 



where T) 



(fc) 



irp(fc) r(fc)i 



Each diagonal block of Eq. 



142] is just the form the 2k + 1 D LL problem of non- 
relativistic fermions in Ref. 

Again we can define the following conserved quantity 



K = 



(D-2)h 











h) 



(37) 

K divides the eigenstates into two sectors flK t 



^ •.<•! in ! ' 



h(l + D 



(43) 



•,3>{im} 



(44) 



respectively. As explained in Sect. IIV[ j represents the 
spin-orbit coupled representation for the SO(D = 2k + 
1) group, and {jm} represents a set of good quantum 
number of the subgroup chain from SO(D — 1) down to 
50(2). 

Similarly as before, the sectors of j ^ j are dis- 
persionless and dispersive with respect to j, respectively. 
The concrete wavefunctions are the same as those in Eq. 
l25l and Eq. [27] by replacing the 3D wavefunction to the 
D-dimensional version of Eq. 1391 The wavefunctions of 



i, ii 



±n r . 3 ,{j m } are § iven explicitly as 

(r) 



±n r ,j,j z 1 



1 

71 



(r) 

l,S,J_,m,Dm}^ 



Y n r ,s,j + ,l,{j m } 



77 

±n r ,j 



V , „ r , SJ+ ,i,{ Jm }( r; ') 



(0 



(45) 



The dispersion relation for the LL branch of <lH „ • r ■ i 
still behaves as E± nr .j^j z y = ±hu^yn r , while that 
for the branch of $!± n j r j j reads as E± 7lr .j {j^ — 



n r + l + %. 



Again each occupied zero energy LL contributes 
to i-fcrmion vacuum charge. If the zeroth LLs 
are fully filled, the vacuum charge is still p D (r) 
2 Sj,{j m } \^0,j,{j m } 



(r)\ 2 , which is expressed as 



? {gr(; + f)(; 



(46) 
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where D — 2k + 1; gi(D) is the degeneracy of the pos- 
itive helicity sector of the fundamental spinor coupling 
to the l-th .D-dimensional spherical harmonics, and its 
expression is the same as dj + given in Eq. 1371 Qd = 
Dw D / 2 /T(D/2 + 1) is the area of D-dimensional unit 
sphere. Eq. [46] can be summed analytically as 



p(r) = 



y/2( 
4 



Wig 



F(D-1, 



D 



(27r)^r(^i) 



< j2 r 

D-2 



as r 



oo. 



(47) 



Similarly, if the /^-dimensional version of the mass 
term inside Eq. [§] is added, every wavefunction with 
the radial quantum number n r hybridizes with its part- 
ner with —n r while keeping all other quantum num- 
ber s the same. T he pair of new eigenvalues becomes 
±y/ (huj) 2 n r + A 2 . Again, the zero-th LL wavefunctions 
single out and remain the same, but their energies are 
shifted to A. For a similar open surface problem to that 
in Sect. MI C| each LL with n r ^ develops a branch of 
gapless surface mode with the upturn (downturn) disper- 
sion with respect to j for n r > (n r < 0), respectively. 
The surface mode from the zeroth LL develops either up- 
turn or downturn dispersions depending on the relative 
sign of the background field coupling and the vacuum 
mass. 



Both of H± commute with the total angular momentum 
operators in the combined representation of s©s' defined 



J i 7 : 



ij;s(Bs' 



Li 



Si- 



(50) 



We choose H + as an example to present the solutions 
of the LL wavefunctions in even dimensions. The K- 
operator is similarly defined as in Eq. [43] as 



K, = 



2S ij L ij + (D - 2)h 






-(2S i jLij 



h) 



, (51) 



where i,j run from 1 to 2k, and Sij and S[j are generators 
in the two fundamental spinor representations given in 
Eqs. [35l and l36l respectively. They again can be divided 
into two sectors of fy + < 1 and <Ir+> 7/ whose eigenvalues of 
K+ are h(l + D — 2) and —HI, respectively. The dis- 
persionless branch of ^P+' 7 can be viewed as LLs, whose 
wavefunctions read 



(f) 



1 
71 



VVj + ,s,Z,{ Jm }(r) 
= F#r»r-l,3-,«',J+l,{jm}(*0 



(52) 



Their spectra are same as before ^±n r ,j,{ i j m } = ihu> 
Please note that the upper and lower components involve 
the s and s' representations, respectively. Similarly, the 
dispersive solutions of "J 77 become 



+Ji 



1 

71 



(53) 



VI. THE LLS OF EVEN DIMENSIONAL 
DIRAC FERMIONS 

The LL problem in the even dimensions with D — 2k is 
more complicated. The SO(2k) group has two irreducible 
fundamental spinor representations s and s'. Each of 
them is with the dimension of 2 k ~ 1 . The construction of 
the SO(2k) generators for the irreducible representations 
are introduced in Sect. IIVI 

(k) ~ 

Now we define i?2fe — o-i where i runs over I to 

2k. Similarly to the 2D case, the counterpart of Eq. [7] 
in the D = 2k dimensions Ho = 2k is reducible into two 



2 blocks as 





±a 2 k - ir 



(fe-l). 



±4fc 







, (48) 



where the repeated index i runs over from 1 to 2k — 1. 
For each one of the reduced Hamiltonians H± , each off- 
diagonal block has only the SO(2k — I) symmetry. Nev- 
ertheless, each of H± is still SO(2k) invariant. If we com- 
bine the two irreducible fundamental spinor representa- 
tions s and s' together, the spin generators are defined 
as 



Si 



A [ i ' 3 J 



(49) 



whose dispersions read -E±n r .j,{j m } = ±fio;yn r + I + -j. 
The solutions to are very similar to Eq. [52] and Eq. 
I53l bv exchanging the irreducible fundamental spinor rep- 
resentation indices s and s'. 

Again if the zeroth branch LLs are filled, the vacuum 
charge p D {f) = \Y,Ui \ 1*0 i U \^)\ 2 is calculated as 



p D (r) 



1 °° 1 r 



1 9l(D)\ -f 



n, 



}e (.VI) 



where D = 2k, and gi{D) are defined similarly as in 
Eq. 2ni It can be summed over analytically as 



p D (r) 



which p D (r) 



4 Vtt/o 



' D r 2 - 



(55) 



D-2 



as r 



VII. CONCLUSION 



In summary, we have generalized the LL problem of 2D 
Dirac fermions to arbitrary higher dimensional flat spaces 
with spherical symmetry. This problem is essentially the 
square root problem of its non-relativistic LL problem 
with spherical symmetry in high dimensions Ref. [l|. The 



8 



zero energy LLs is a branch of i-fermion modes. On the 
open boundary, each LL contributes one branch of helical 
surface modes. This series of LL problems can be viewed 
as the generalization of parity anomaly in 2D to arbitrary 
dimensions in a spherical way. An open question is that 
how to experimentally realize the case of the 3D systems. 

Note added Near the completion of this manuscript, 
we became aware of that Eq. [5] plus Eq. [9] with a special 
relation between the background field coupling constant 



and the mass term was studied in Ref. (4lU42| under 
the name of Dirac oscillator. We have studied a more 
general form from a different perspective and identified 
their relation to the LLs. 

This work was supported in part by the NBRPC (973 
program) 2011CBA00300 (2011CBA00302) and Grant 
No. ARO-W911NF0810291. K. I. was supported by 
DOE-FG03-97ER40546. 
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